Hybrid zone formation and contrasting outcomes of secondary contact over transects in common toads.

Much progress in speciation research stems from documenting patterns of morphological and genetic variation in hybrid zones. Contrasting patterns of marker introgression in different sections of the contact can provide valuable insights on the relative importance of various evolutionary mechanisms maintaining species differences in the face of hybridization and gene flow and on hybrid zone temporal and spatial dynamics. We studied species interactions in the common toads Bufo bufo and B. spinosus in France and northwestern Italy using morphological and molecular data from the mitochondrial and nuclear genomes in an extensive survey, including two independent transects west and east of the Alps. At both, we found sharp, coincident and concordant nuclear genetic transitions. However, morphological clines were wider or absent and mtDNA introgression was asymmetric. We discuss alternative, nonexclusive hypotheses about evolutionary processes generating these patterns, including drift, selection, long-distance dispersal and spatial shifts in hybrid zone location and structure. The distribution of intraspecific mtDNA lineages supports a scenario in which B. bufo held a local refugium during the last glacial maximum. Present-day genetic profiles are best explained by an advance of B. spinosus from a nearby Iberian refugium, largely superseding the local B. bufo population, followed by an advance of B. bufo from the Balkans, with prongs north and south of the Alps, driving B. spinosus southwards. A pendulum moving hybrid zone, first northwards and then southwards, explains the wide areas of introgression at either side of the current position of the contact zones.

Molecular phylogenetics and historical biogeography of the west-palearctic common toads (Bufo bufo species complex).

In most pan-Eurasiatic species complexes, two phenomena have been traditionally considered key processes of their cladogenesis and biogeography. First, it is hypothesized that the origin and development of the Central Asian Deserts generated a biogeographic barrier that fragmented past continuous distributions in Eastern and Western domains. Second, Pleistocene glaciations have been proposed as the main process driving the regional diversification within each of these domains. The European common toad and its closest relatives provide an interesting opportunity to examine the relative contributions of these paleogeographic and paleoclimatic events to the phylogeny and biogeography of a widespread Eurasiatic group. We investigate this issue by applying a multiproxy approach combining information from molecular phylogenies, a multiple correspondence analysis of allozyme data and species distribution models. Our study includes 304 specimens from 164 populations, covering most of the distributional range of the Bufo bufo species complex in the Western Palearctic. The phylogenies (ML and Bayesian analyses) were based on a total of 1988 bp of mitochondrial DNA encompassing three genes (tRNAval, 16S and ND1). A dataset with 173 species of the family Bufonidae was assembled to estimate the separation of the two pan-Eurasiatic species complexes of Bufo and to date the main biogeographic events within the Bufo bufo species complex. The allozyme study included sixteen protein systems, corresponding to 21 presumptive loci. Finally, the distribution models were based on maximum entropy. Our distribution models show that Eastern and Western species complexes are greatly isolated by the Central Asian Deserts, and our dating estimates place this divergence during the Middle Miocene, a moment in which different sources of evidence document a major upturn of the aridification rate of Central Asia. This climate-driven process likely separated the Eastern and Western species. At the level of the Western Palearctic, our dating estimates place most of the deepest phylogenetic structure before the Pleistocene, indicating that Pleistocene glaciations did not have a major role in splitting the major lineages. At a shallow level, the glacial dynamics contributed unevenly to the genetic structuring of populations, with a strong influence in the European-Caucasian populations, and a more relaxed effect in the Iberian populations.

Multilocus species tree analyses resolve the radiation of the widespread Bufo bufo species group (Anura, Bufonidae).

New analytical methods are improving our ability to reconstruct robust species trees from multilocus datasets, despite difficulties in phylogenetic reconstruction associated with recent, rapid divergence, incomplete lineage sorting and/or introgression. In this study, we applied these methods to resolve the radiation of toads in the Bufo bufo (Anura, Bufonidae) species group, ranging from the Iberian Peninsula and North Africa to Siberia, based on sequences from two mitochondrial and four nuclear DNA regions (3490 base pairs). We obtained a fully-resolved topology, with the recently described Bufo eichwaldi from the Talysh Mountains in south Azerbaijan and Iran as the sister taxon to a clade including: (1) north African, Iberian, and most French populations, referred herein to Bufo spinosus based on the implied inclusion of populations from its type locality and (2) a second clade, sister to B. spinosus, including two sister subclades: one with all samples of Bufo verrucosissimus from the Caucasus and another one with samples of B. bufo from northern France to Russia, including the Apennine and Balkan peninsulas and most of Anatolia. Coalescent-based estimations of time to most recent common ancestors for each species and selected subclades allowed historical reconstruction of the diversification of the species group in the context of Mediterranean paleogeography and indicated a long evolutionary history in this region. Finally, we used our data to delimit the ranges of the four species, particularly the more widespread and historically confused B. spinosus and B. bufo, and identify potential contact zones, some of which show striking parallels with other co-distributed species.

Molecular phylogenetic relationship of toads distributed in the Far East and Europe inferred from the nucleotide sequences of mitochondrial DNA genes.

The toads of the Bufo bufo species group are widely distributed in the Eurasian continent and Japanese Archipelago. In this study, we analyzed the mtDNA gene sequences of this species group and estimated the divergence time to clarify the evolutionary relationships and biogeography of toads distributed in the Far East and Europe. The phylogenetic tree indicated that this group produced Bufo bufo in Europe, whereas it produced B. japonicus in the Far East. B. japonicus was divided into three major clades corresponding to a group consisting of B. j. gargarizans in China, B. j. bankorensis in Taiwan, and B. j. miyakonis on Miyako Isl. and eastern and western groups of Japanese B. j. japonicus subspecies group. The eastern and western groups were divided into several subclades which tended to reflect the region-specific geographic distribution of all localities except B. j. japonicus from Hakodate. The estimated branching times of these clades suggest that geological events may have influenced the divergence of the toads distributed in the Far East and Europe.